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Abstract 

Motivated by the recent LHC hints of a Higgs boson around 125 GeV, we assume a SM-hke Higgs 
with the mass 123-127 GeV and study its implication in low energy SUSY by comparing the MSSM 
and NMSSM. We consider various experimental constraints at 2a level (including the muon g — 2 
" and the dark matter relic density) and perform a comprehensive scan over the parameter space of 

each model. Then in the parameter space which is allowed by current experimental constraints and 
also predicts a SM-like Higgs in 123-127 GeV, we examine the properties of the sensitive parameters 
(like the top squark mass and the trilinear coupling At) and calculate the rates of the di-photon 
signal and the VV* {V = W, Z) signals at the LHC. Our typical findings are: (i) In the MSSM the 
top squark and At must be large and thus incur some fine-tuning, which can be much ameliorated 
in the NMSSM; (ii) In the MSSM a light stau is needed to enhance the di-photon rate of the 
SM-like Higgs to exceed its SM prediction, while in the NMSSM the di-photon rate can be readily 
enhanced in several ways; (iii) In the MSSM the signal rates of — )• /i — )• VV* at the LHC are 
never enhanced compared with their SM predictions, while in the NMSSM they may get enhanced 
significantly; (iv) A large part of the parameter space so far survived will be soon covered by the 
expected XENON100(2012) sensitivity (especially for the NMSSM). 

PACS numbers: 
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I. INTRODUCTION 



As the only missing particle in the Standard Model (SM), the Higgs boson is being 
intensively hunted at the LHC Recently, both the ATLAS and CMS experiments have 
revealed hints of a Higgs particle around 125GeV 

HQ 

. While such a Higgs mass can be 
well accommodated in the SM and the reported signal rates in several channels are also 
in agreement with the SM expectations after taking into account the large experimental 
uncertainty 3|, \^ (albeit the central value of the observed di-photon rate is somewhat above 
the SM prediction), the low energy supersymmetry (SUSY) seems to be a better framework 
to account for such a Higgs. In low energy SUSY the SM-like Higgs mass is theoretically 
restricted in a narrow range and its di-photon rate at the LHC may exceed the SM prediction 
5|-[7|, both of which are welcomed by the LHC results. 

However, as the most popular low energy SUSY model, the minimal supersymmetric stan- 
dard model (MSSM) js,!^ may have some tension to accommodate such a 125 GeV Higgs. 
As is well known, in the MSSM the SM-like Higgs mass is upper bounded by Mz cos 2/3 at 
tree- level and to get a Higgs around 125GeV we need sizable top/stop loop contributions, 
which depend quartically on the top quark mass and logarithmically on the stop masses 
[lo| . This will, on the one hand, impose rather tight constraint on the MSSM, and on the 
other hand, incur some fine-tuning Such a problem can be alleviated in the so-called 
next-to-minimal supersymmetric standard model (NMSSM) 12, ll3|, which is the simplest 
singlet extension of the MSSM with a scale invariant superpotential. In the NMSSM, due to 
the introduction of some new couplings in the superpotential, the SM-like Higgs mass gets 
additional contribution at tree-level and also may be further pushed up by the mixing effect 
in diagonalizing the mass matrix of the CP-even Higgs fields 14|, [l5|. As a result, a SM-like 
Higgs around 125GeV does not entail large loop contributions, which may thus ameliorate 



the fine-tuning problem 



16|. 



In this work, motivated by the recent LHC results, we assume a SM-like Higgs boson 
in 123 - 127 GeV and study its imphcation in the MSSM and NMSSM. Different from 



recent studies in this direction 
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18| . we scan the model parameters by con- 



sidering various experimental constraints and perform a comparative study for the MSSM 
and NMSSM. We will investigate the features of the allowed parameter space in each model 
and particularly we pay more attention to the space of the NMSSM which may be dis- 
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tinct from the MSSM. Noting the LHC experiments utihze the channels pp h ^ 
pp ^ h ^ ZZ* —7- 4Z and pp ^ h ^ WW* — ?■ 21 2v in searching for the Higgs boson 
we also study their normalized production rates defined as 



^77 = (^susy{pp -^h^ 1i)/(^sm{pp ^ h ^ j-f) 

^ Clfl, X TUhsM)/TUh). (1) 

Rvv = (ysusvipp VV*)/asM{pp ^ h VV*) 

= ClggClyy X TUhsM)/rtotih), (2) 

where V = W, Z, and Chgg, Ch'y-y and Chvv are respectively the rescaled couplings of the 
Higgs to gluons, photons and weak gauge bosons by their SM values. We are more interested 



in the case with R^^ > 1 because it is favored by current ATLAS and CMS results [liJal. As 
will be shown below, this case usually predicts a slepton or a chargino lighter than 250 GeV. 

This paper is organized as follows. In Sec. II, we recapitulate the characters of the 
Higgs mass in the models for better understanding our numerical results. In Sec. Ill, we 
perform a comprehensive scan over the parameter space of each model by imposing current 
experimental constraints and also by requiring a SM-like Higgs boson in 123 — 127GeV. 
Then we scrutinize the properties of the surviving parameter space. Finally, we draw our 
conclusions in Sec. IV. 

II. THE SM-LIKE HIGGS MASS IN THE MSSM AND NMSSM 

In the MSSM, the Higgs sector consists of two doublet fileds and H^, which after the 
electroweak symmetry breaking, result in five physical Higgs bosons: two CP-even scalars h 
and H, one CP-odd pseudoscalar A and a pair of charged scalars if^[8|. Traditionally, such 
a Higgs sector is described by the ratio of the Higgs vacuum expectation values, tan (3 = 
and the mass of the pseudoscalar m^. In most of the MSSM parameter space, the lightest 
Higgs boson h has largest coupling to vector bosons (i.e. the so-called SM-like Higgs boson), 



and for moderate tan /3 and large its mass is given by 18| 



mi - M| cos' 2/j + — ^ In ^ + ^ , \ Jo [ ^ - TTTTT? > (3) 



where v = 174 GeV, Ms = ^Jmi^mi^ with m^^ and m^^ being the stop masses, Xt = 
At — 11 cot P with At denoting the trilinear Higgs-stop coupling and being the Higgsino 



mass parameter. Obviously, the larger tan /3 or Ms is, the heavier h becomes, and for given 



M5, m/j reaches its maximum when XtjMg = VQ, which corresponds to the so-called m^"^ 



scenario. 



About Eq.(l3]), three points should be noted [18|. The first is this equation is only valid 
for small splitting between m^-^ and m^^. In case of large splitting, generally Xt/Ms > V6 
is needed to maximize mh. The second is ml in Eq.(|3]) is symmetric with respect to the 
sign of Xf. This behavior will be spoiled once higher order corrections are considered, and 
usually a larger ruh is achieved for positive AtM^ with M3 being gluino soft breaking mass. 
And the last is in Eq.OH]), we do not include the contributions from the sbottom and slepton 
sectors. Such contributions are negative and become significant only for large tan /3. 

Compared with the MSSM, the Higgs sector in the NMSSM is rather complex, which can 
be seen from its superpotential and the corresponding soft-breaking terms given by [l2| 

W^NMSSM = Wf + XHu ■ HdS + ^^S^ (4) 

P.NMSSM ^ ^2 |^^|2 ^ ^2|^^|2 + ^2 |^|2 ^ ^A.XSH^ ■ H, + ^kS' + h.c). (5) 

Here Wp is the superpotential of the MSSM without the /i term, the dimensionless param- 
eters A and k are the coefficients of the Higgs self couplings, and m„, m^, m^. Ax and 
are the soft-breaking parameters. 

After the electroweak symmetry breaking, the three soft breaking masses squared for 
Hu, Hd and S can be expressed in terms of their VEVs (i.e. t>„, Vd and s) through the 
minimization conditions of the scalar potential. So in contrast to the MSSM where there 
are only two parameters in the Higgs sector, the Higgs sector of the NMSSM is described 
by six parameters [l^ : 

A, K, = — , A^, tan/3 = — , /i = As. (6) 

sm 2p Vd 

The Higgs fields can be written in the following form: 

H,= \ J'^ I, H,= \ I, Hs = s + ^iS, + ^P2), (7) 















Sl+iPl 
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where Hi = cos f3Hu — e sin I3H^, H2 = sin (3Hu + e cos f3H^ with £12 = £21 = —1 and 
= £^22 = 0, and are Goldstone bosons and v = \/vl + f^. In the CP-conserving 
NMSSM, the fields 5*1, 5*2 and 5*3 mix to form three physical CP-even Higgs bosons, and Pi 



and P2 mix to form two physical CP-odd Higgs bosons. Obviously, the field H2 corresponds 
to the SM Higgs field, and the scalar h with largest S2 component is called the SM-like Higgs 
boson. 

Under the basis (5*1, S2, S-^), the elements of the mass matrix for Si fields at tree level 



are given by 



15 



19| 



11 



MI2 



Ml, 
Ml, 



Mi + (m|-AV)sin^2/3 
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+ + ^) - ^^'^^^ 2/3, 



(8) 
(9) 

(10) 

(11) 
(12) 

4 . - . ^ . . ^ . ^ . 2 ' ^^^^ 

where A^22 is nothing but ml at tree level without considering the mixing among Si, and its 
second term A^f^sin^2/3 originates from the coupling XH^ ■ HiiS in the superpotential. For 
such a complex matrix, it is useful to consider two scenarios for understanding the results: 

• Scenario I: A, k — )■ and /i is fixed. In this limit, since Ail,,J^23 — O5 singlet 
field 5*3 is decoupled from the doublet fields, and the MSSM mass matrix is recovered 
for the (5*1, 5*2) system. This scenario indicates that, even for moderate A and k, rrih 
should change little from its MSSM prediction. So in order to show the difference of 
the two models in predicting m^, we are more interested in large A case. Especially 
we will mainly discuss A > Mz/v ^ 0.53 case, where the tree level contributions to 
m|, i.e. A^225 maximized for moderate values of tan/3 rather than by large values 
of tan /3 as in the MSSM. 

• Scenario II: Mf^ > A^22 > -^12 and {Mji - Ml,) > Mj,, which can be easily 
realized for a large M\. In this limit. Si is decoupled from the {S2,S3) system and the 
properties of rrih can be qualitatively understood by the 2x2 matrix |l5| 



'22 
Ml, 



Ml, - 



(14) 



where 5^ denotes the radiative corrections to rrih with its form given by the last two 
terms of Eq.(l3]), and represents the potentially important effect of (5*1, 5*3) mixing on 



Mig. This matrix indicates that the ( 5*2, 5*3) mixing can push rrih up once Mf^ > MI2, 
and such effect is maximized for Mf-^ shghtly larger than and at the same time 



M^2 = -^23 shghtly below Mf^M|2 (larger M12 will destabilize the vacuum) 15|. 
Obviously, in this push-up case, h is the next-to-lightest CP-even Higgs boson and 
the larger the (>S'2,5'3) mixing is, the heavier h becomes. Alternatively, the mixing 
can pull m/j down on the condition of Mf-^ < MI2, which occurs for large Kfi (for 
M^sin2/3/yU ~ 2, see discussion below) as indicated by the expression of A^ls and the 
positiveness of A^. Here we remind that, due to the extra contribution A^f ^ sin^ 2/3 to 
m\ at tree level, rrih in the pull-down case may still be larger than its MSSM prediction 
for a certain 5^. 

Since our results presented below is approximately described by scenario II, we now 
estimate the features of its favored region to predict ^ 125GeV. First, since 
Mfi ~ O(1002)GeV2, Ml^ > M^^ implies that Ma > O(300)GeV. Numerically, we 
find Ma > C(300)GeV for the push-up case, and Ma > C(500)GeV for the pull-down 
case (see Fig|7]). Second, M^2 = -^23 niust be relatively small, which implies that 
MAsm2P/ fi ~ 2 for A > 0.53. This can be understood as follows. In the push-up 
case, since MI2 < ^11 ~ O{100'^)GeV^, the condition A^23 < Mf^M^^ (for vacuum 
stability) has limited the size of A^23- While in the pull-down case, a very large A^23 
will suppress greatly rrih to make it difficult to reach 125GeV, and this in return limits 
the size of A^23- Given that fi > 100 GeV as required by the LEP bound on chargino 
mass and that a larger fi is favored for the pull-down scenario, one can infer that 
the value of sin 2/3//i should be around 2 after considering that the third term 
in AI23 is less important. Numerically speaking, we find {2Xfiv)\ < 0.2 and 

1.4 < Ma sin 2/3 //i < 2 (see Fig]?]). Lastly, light stops may be possible in the NMSSM 
with large A to predict rrih — 125GeV. To see this, we consider the parameters A = 0.7 
and tan/3 = 1.5, and we get 5^/125^ ~ 5% without considering the mixing effect to 
predict rrih — 125GeV. This is in sharp contrast with 5^/125^ ~ 55% in MSSM for 
tan (3 = 5. 



In this work we use the package NMSSMTools 20| to calculate the Higgs masses and 



mixings, which includes the dominant one-loop and leading logarithmic two-loop corrections 



to A^^. We checked our MSSM results of nih by using the code FeynHiggs[2l| and found the 
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results given by NMSSMTools and FeynHiggs are in good agreement ( for ruh ~ 125GeV 
they agree within O.bGeV for same MSSM parameters). 



III. NUMERICAL RESULTS AND DISCUSSIONS 

In this work, we scan the parameters of the models and investigate the samples that 
predict 123GeV < rrih < 127GeV and at the same time survive the following constraints 



20( |: (1) The constraint from the LHC search channel pp H ^ 2t for non-standard 
Higgs boson. (2) The limits from the LEP and the Tevatron on the masses of sparticles 
as well as on the neutralino pair productions. (3) The constraints from B-physics, such as 
B — )■ Xg'y, the latest experimental result of Bg — ?■ fi^fJ'^ , Bd — ?■ Xg/i+yU^, B^ r+z/ and the 
mass differences AM^ and AM^. (4) The indirect constraints from the electroweak precision 



observables such as siv?6l^j^, pi and My[r, and their combinations ei{i = 1,2,3) 



3- We 

require to be compatible with the LEP/SLD data at 95% confidence level. We also require 
the SUSY prediction of the observable Rb (X{Z — )• bb)/T{Z — )■ hadrons)) to be within the 2a 
range of its experimental value 23|. (5) The constraints from the muon anomalous magnetic 



moment: al^'P - a^^'' = (25.5 ± 8.2) x 10^^° We require the SUSY effects to explain the 



discrepancy at 2a level. (6) The dark matter constraints from WMAP relic density (0.1053 
< nh'^ < 0.1193) I25I and the direct search result from XENONIOO experiment (at 90% 



CL.) 26|. (7) For the NMSSM, we also require the absence of a landau singularity below 



the GUT scale, which implies A < 0.7 for small k and k < 0.5 for A > 0.53 at weak scale. In 



our calculation, we fix rrit = 172.9 GeV and frs = 0.02 [27| (/ts denotes the strange quark 
fraction in the proton mass), and use the package NMSSMTools to implement most of the 
constraints and to calculate the observables we are interested in. 

In our scan, we note that the soft parameters in the slepton sector can only affect signif- 
icantly the muon anomalous magnetic moment a^, which will in return limit the important 
parameter tan/3, so we assume them a common value m|- and treat it as a free parameter. 
For the soft parameters in the first two generation squark sector, due to their little effects 
on the properties of the Higgs bosons, we fix them to be 1 TeV. As for the gaugino masses, 
we assume the grand unification relation, 3Mi/5ai = M2/a2 = M^/a^ with being the 
fine structure constants of the different gauge groups, and treat Mi as a free parameter. In 
order to reduce free parameters, we also assume the unimportant parameters and Af, 
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to satisfy M^g = M^c/g and Ab = At. 



A. Implication of rrih — 125GeV in generic SUSY 

In order to study the implication of rrih — 125GeV in generic SUSY, we relax the soft mass 
parameters to 5 TeV and perform an extensive random scan over the following parameter 
regions: 

1 < tan /3 < 60, 90 GeV < < 1 TeV, 100 GeV < < 2 TeV, 
100 GeV < Mq3, Mc/3 < 5 TeV, \At\ < 5 TeV, 

100 GeV <mi<l TeV, 50 GeV < Mi < 500 GeV, (15) 
for the MSSM, and 

< A < 0.2, 0<K< 0.7, 90 GeV < Ma < I TeV, \A^\ < 1 TeV, 
100 GeV < Mq3, M(/3 < 5 TeV, \At\ < 5 TeV, 

1 < tan /3 < 60, 100 GeV < fi, mj < 1 TeV, 50 GeV < Mi < 500 GeV, (16) 

for the NMSSM. In our scan, we only keep the samples satisfying the requirements listed 
in the text (including 123GeV < rrih < 127GeV). To show the differences between the 
MSSM and the NMSSM, we also perform a scan similar to Eq. (fT6|) except that we require 
A > rriz/v ~ 0.53. In FigJH we show the correlation of the lighter top-squark mass (?tijJ with 
the ratio Xt/Ms {Ms = ^/rr^Jri^) for the surviving samples in the MSSM and NMSSM. 
As expected from Eq.(l3]), in order to predict rrih — 125GeV in the MSSM, a large Xt is 
needed for a moderate light ti, and with ti becoming heavy, the ratio Xt/Ms decreases, 
but is unlikely to vanish for rrii^ < 5TeV. These features are maintained for NMSSM with 
A < 0.2 (see the middle panel) but changed for NMSSM with a large A, where Xt can 
possibly vanish even for rrti^ ~ ITeV. FigH] also shows that a ti as light as 200GeV is still 
able to give the required rrih. But in this case Xt is large [Xtj ^mi^rrii^ > a/6), which leads 
to a large mass splitting between two stops (mj^ ^ "^ti)- Note that a ii as light as 200GeV 
does not contradict the recent SUSY search result of the LHC 28|. 

Since in heavy SUSY the radiative correction 6"^ is usually very large, rrih — 125GeV is 
unlikely to impose tight constraints on other parameters of the models. Considering heavy 
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^•"^3 X,/Ms X,/Mg 

FIG. 1: The scatter plots of the samples in the MSSM and NMSSM satisfying all the requirements 
(1-7) listed in the text (including 123GeV < < 127GeV), projected in the plane of m^^ versus 
Xt/Ms with Ms = y/f^ti "^te At — /icot j3. 



SUSY is disfavored by naturalness, we in the following concentrate on the implication of 
nih ^ 125GeV in sub-TeV SUSY. 



B. Implication of mn ~ 125GeV in sub-TeV SUSY 

In this section, we study the implication of rrih — 125GeV in low energy MSSM and 
NMSSM. In order to illustrate the new features of the NMSSM, we only consider the case 
with A > 0.53. Our scans over the parameter spaces are quite similar to those in Eq. (fT5|) 
and Eq.(lT6l) except that we narrow the ranges of Mgg, Mjj^ and At as follows: 

100 GeV < (Mq3, M[/3) < 1 TeV, | A| < 3 TeV. (17) 

In Fig|2] we project the surviving samples of the models in the plane of m^^ versus At, 
showing the results with i?^^ < 1 and i?^^ > 1 separately. As we analyzed in Sec. II, the 
SM-like Higgs in the NMSSM may be either the lightest Higgs boson (corresponding to the 
pull-down case) or the next-to-lightest Higgs boson (the push-up case). In the figure we 
distinguished these two cases. We note that among the surviving samples the number of the 
pull-down case is about twice the push-up case. 

Figj2]shows that, in order to get mh — 125GeV in the MSSM, mi^ and \At\ must be larger 
than about 300GeV and 1.5TeV respectively, and the bounds are pushed up to 600GeV and 
l.STeV respectively for R^^ > 1. While in the NMSSM, a ii as hght as about lOOGeV (in 
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FIG. 2: Same as Fig.l, except that the samples are based on the narrowed scan ranges of the soft 
masses shown in Eq. (ll7p for both models and the requirement A > 0.53 for the NMSSM. Here 
the samples are projected in the plane of m^^ versus Af. The upper (lower) panels correspond to 
R-y-y < 1 {R.yy > l) wlth R^^ = (T su SY {pp k ^ n) / s M {pp k ^ . For the NMSSM 
results, the circles (green) denote the case of the lightest Higgs boson being the SM-like Higgs (the 
so-called pull-down case), and the times (red) denotes the case of the next-to-lightest Higgs boson 
being the SM-like Higgs (the so-called push-up case). 

either the pull-down or push-up case) is still able to predict ^ 125GeV, and even if one 
require R^^ > 1, a ti as light as about 200GeV is allowed. The fact that the NMSSM allows 
a lighter ti than the MSSM indicates that the NMSSM is more natural than the MSSM in 
light of the LHC results. 

Since a light ti may significantly change the effective couplings Chgg and Ch-^^, we present 
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TABLE I: The ranges of the rescaled couphngs Chgg = CfUSY/^SM (Jh^^ = ^SUSY/^SM 
predicted by the surviving samples of the two models. The region of C^ggiR-y-y > 1) is obtained by 
only considering the samples with R^^ > 1. 
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FIG. 3: Same as Figl2l but showing the dependence of the di-photon signal rate R-yy on the effective 
hbb coupling = C7SUSY/C7SM_ 



in Table I their predicted ranges for the surviving samples. This table shows that Chgg is 
always reduced, and for m^^ ~ lOOGeV in the NMSSM, the reduction factor may reach 70%. 
While Ch-y-y exhibits quite strange behaviors: it is enhanced in the MSSM, but may be either 
enhanced or suppressed in the NMSSM. This is because, unlike Chgg which is affected only 
by squark loops, Ch^yy gets new physics contributions from loops mediated by charged Higgs 
boson, charginos, sleptons and also squarks, and there exists cancelation among different 
loops. As will be shown below, the current experiments can not rule out light sparticles like 
fi and chargino. Although the contributions of these particles to Ch-yy are far smaller than 
the W loop contribution, they may still alter the coupling significantly. 

Since the di-photon signal is the most important discovery channel for the Higgs boson 
around 125GeV, it is useful to study its rate carefully. From Eq.([T]) one can learn that the 
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rate is affected by Chgg and Ch-y^ discussed above, and also by the total width of h (or more 
basically by the hhh coupling since hh is the dominant decay of h). The importance of hhh 
coupling on the di-photon rate was recently emphasized in 6[]. Here we restrict our study 
to the nih — 125GeV case. In Fig|31 we show the dependence of the di-photon signal rate 
R^^ on the effective hhb coupling (including the potentially large SUSY corrections {g, Q) 
normalized by its SM value. This figure indicates that although the rate is suppressed for 
most of the surviving samples in both models, there still exist some samples with enhanced 
rate, especially the NMSSM is more likely to push up the rate than the MSSM. This feature 
can be understood as follows. In SUSY, i?^^ > 1 requires approximately the combination 
ChggCh-yy/Chbb cxcccd 1. For the MSSM, given Chgg < 1 and C^bb > 1 for nearly all 
the cases, this condition is not easy to satisfy. While in the NMSSM, due to the mixing 
between the doublet field 5*2 and the singlet field 5*3, Chbh < 1 is possible once the singlet 
component in h is significant, which is helpful to enhance the combination. In fact, we 
analyzed carefully the R-y^ > 1 cases and found they are characterized by Chgg, Chbi — 1 and 
Ch-y-y > 1 in the MSSM, and by C^bi < 1 in the NMSSM. In other words, it is the enhanced 
/177 coupling (reduced total width of h) that mainly push up the di-photon rate in the 
MSSM (NMSSM) to exceed its SM prediction. Figj3]also indicates that the pull-down case 
in the NMSSM is less effective in reducing Chhi and thus can hardly enhance the di-photon 
rate. This is because in the push-up case, both M^^ and in Eq. ffMl) are moderate and 
often comparable, which are helpful to enhance the (5*2, 6*3) mixing. Finally, we note that 
in some rare cases of the NMSSM the ratio R^^ may be very small even for Chbi < 1- This 
is because there exists very light Higgs boson so that h decays dominantly into them. 

In order to further clarify the reason for the enhancement of the di-photon rate in the 
MSSM, we scrutinize carefully the parameters of the model and find that the samples with 
R^^ > 1 correspond to the case with a large fi tan 13 and m^^ < 200GeV, which is illustrated 
in FigJH This means that the stau loop plays an imp ortant role in enhancing C^^. We note 
that similar conclusion was recently achieved in 18|, but in that work the authors did not 
consider the tight experimental constraints. From FigJH we also note that R^^ < 0.95 is 
predicted in the MSSM with m^^ > 250GeV. So future precise measurement of R^^ and m^^ 
may be utilized to verify the correctness of the MSSM. 

Considering the process pp ^ h ^ VV* {V = W, Z) is another important Higgs search 
channel, we in FiglS] show the signal rate versus the hVV coupling. This figure shows that 
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FIG. 4: Same as Figl^l but only for the MSSM, projected in the planes of tan/5 versus and R. 
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in the MSSM, h is highly SM-like, while in the NMSSM, the singlet component in h may be 
sizable, especially in the push-up case, so that Chvv is reduced significantly. The signal rate 
Rvv also behaves differently in the two models. In the MSSM, because Chgg < 1 and in most 
cases Cfii,i > 1, Rw is always less than 1 (for i?^^ > 1 it varies between 0.7 and 0.95). In 
the NMSSM, however, Rw may exceed 1 and in this case we find Rvv — -R77. The reason 
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FIG. 6: Same as Fig.2, but projected in the planes of /x versus tan /3 and versus k respectively. 

for such a correlation is the two quantities have the same origin for their enhancement, i.e. 
the suppression of the hbb due to the presence of the singlet component in h. 

Next we investigate the favored parameter space of the NMSSM to predict rrih — 125GeV. 
As introduced in Sec. II, besides the soft parameters in the stop sector, the sensitive pa- 
rameters include tan /3, jjj^ Kj clS well as Myi. In Figj6l we project the surviving samples in 
the planes of /i versus tan /3 and versus k. This figure shows three distinctive characters for 
the allowed parameters. The first is that tan /3 must be moderate, below 4 and 9 for the 
pull-down and the push-up case, respectively. Two reasons can account for it. One is that 
in the NMSSM with large A, the precision electroweak data, i.e. the constraint (4), strongly 



disfavor a large tan/3 30|]. The other reason is that, as far as A > 0.53 is concerned, a 
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FIG. 7: Same as FigEl but showing the correlation between and ^/ sin 2/3. The dashed Une 
denotes the relation M^sin2/3//i = 2. 

moderate tan/3 is welcomed to enhance the tree level value of ml (i.e. that even 

without heavy stops, rrih can still reach 125GeV. Moreover, since the (5*2, 5*3) mixing is to 
reduce the value of Mf^ in Eq. lJT^ in the pull-down case, a larger (or equivalently a 
smaller tan/3) is favored by the Higgs mass. The second character is that k/x in the push-up 
case is usually much smaller than that in the pull-down case. This is because, as we intro- 
duced in Sec. II, a large K/i is needed by the pull-down case to enhance MI2 in Eq. (fT^ . The 
third character is obtained by comparing the parameter regions in the upper panels with 
those in the lower panels, which shows that R^^ > 1 puts a lower bound on k, i.e. k > 0.1. 
The underlying reason is that for k < 0.1, the dark matter will be light and singlino-like. 



and to get its currently measured relic densi 
universe by exchanging a light Higgs boson 



2, the dark matter must annihilate in the early 



3l|. In this case, h mainly decays into the light 



bosons, which in return will suppress the di-photon rate. 

In FigJT] we show the correlation of Ma with /x/ sin 2/3. This figure indicates that 
Ma ^ SOOGeV for the push-up case and Ma ^ SOOGeV for the pull-down case, which 
is in agreement with our expectation. In fact, we checked each surviving sample and found 
it satisfies the condition: M| > MI2 > -^12 (^1 ~ -^33) ^ -^i3' ^^e samples 
can be well described by scenario II. We also checked that the mixing of the field 5*1 with 
82/83 is small and Ma is approximately the heaviest CP-even Higgs boson mass. FigJT] also 
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FIG. 8: Same as FigUl but projected on the plane of the charged Higgs boson mass and tan /3. 



indicates that the relation M4sin2/3//i = 2 is maintained quite well in the push-up case, 
but is moderately spoiled in the pull-down case. The reason is, as we introduced in Sec. 
II, the requirement that A^23 should be moderately small actually implies Ca ~ with 
Ca = 1 — (Mii|HLM)2 _ ^ sin 2/3. In the push-up case the third term in Ca is not important, 
while in the pull case, although it is several times smaller than the second term, it may not 
be negligible. We checked our results and found \Ca\ ^ 0.2 and 1.4 < M4sin2/3//i < 2 for 
all the surviving samples. 

About the NMSSM with A > 0.53, three points should be noted. The first is, from the 
results presented in Figj6l one may find the presence of a smuon and/or a chargino lighter 
than 250GeV. This is because the surviving samples are characterized either by tan/3 < 4 
or /i < 250GeV or by both in the NMSSM (see Fig|6]). Then to explain the discrepancy of 
muon anomalous magnetic moment < 250GeV is needed for a low tan /3, and m<^± ^ /i 
implies m^± < 250GeV. We numerically checked the validity of this conclusion. The second 
point is, although M4sin2/3//i ~ 2 may be regarded as a new source of fine tuning in the 
NMSSM, it is rather predictive to get the value of Ma once /i and tan /3 are experimental 
determined. Finally, we note the favored region for fi and tan /3 shown in Figl6] does not 
overlap with that in FiglH This may be used to discriminate the models. 

Finally, we briefly describe other implications of rrih — 125GeV in the SUSY models. In 
FiglHlwe project the surviving samples on the plane of tan/3 versus mH+ with denoting 
the charged Higgs boson. This figure shows that must be heavier than about 200GeV 
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FIG. 9: Same as Figl21 but exhibiting the spin-independent x-nucleon scattering cross section as 
a function of the dark matter mass. 

in the MSSM. This bound is much higher than the corresponding LEP bound, which is 
about 80GeV. For the NMSSM with a large A, the bound can be further pushed up to 
about 300GeV. This figure also indicates that in the MSSM, tan /3 may reach 35 for 771^+ = 
400GeV. Then based on the MC simulation by the ATLAS collaboration 32|, one may 
expect that the charged Higgs may be observable from the process pp — )■ tH~ — )■ bWrUr at 
the early stage of the LHC However, this may be impossible. The reason is, for relatively 
light and large tan/3, /i must be large to satisfy the constraint from dark matter direct 
detection experiments such as XENONIOO. This will greatly suppress the ibH'^ coupling 



331]. For the NMSSM, the hope to observe is also dim because tan/3 is small. 

In Figini we show the spin-independent elastic scattering between dark matter and nu- 



cleon. We use the formula presented in the Appendix of 3J] to calculate the scattering rate. 
As expected, the XENONIOO (2012) data to be released in near future will further exclude 
some samples, especially the pull-down case of the NMSSM will be strongly disfavored if 
XENONIOO (2012) fails to find any evidence of dark matter (assuming the grand unification 
relation of the gaugino mass). From the left of Figj9]one can learn that for the samples with 
R^^ > 1 in the MSSM, the scattering rate is small, usually at least one order below than 
the sensitivity of the XENONIOO (2012). 
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IV. CONCLUSION 



Motivated by the recent LHC hints of a Higgs boson around 125 GeV, we assume a 
SM-hke Higgs with the mass 123-127 GeV and study its imphcation in low energy SUSY 
by comparing the MSSM and NMSSM. Under various experimental constraints at la level 
(including the muon g — 2 and the dark matter relic density), we scanned over the param- 
eter space of each model. Then in the parameter space allowed by current experimental 
constraints and also predicting a SM-like Higgs in 123-127 GeV, we examined the properties 
of the sensitive parameters and calculated the rates of the di-photon signal and the W* 
(y = W, Z) signals at the LHC. Among our various findings the typical ones are: (i) In 
the MSSM the top squark and At must be large and thus incur some fine-tuning, which 
can be much ameliorated in the NMSSM; (ii) In the MSSM a light f is needed to make 
the di-photon rate of the SM-like Higgs exceed its SM prediction, while the NMSSM has 
more ways in doing this; (iii) In the MSSM the signal rates of pp h ^ VV* at the LHC 
are never enhanced compared with their SM predictions, while in the NMSSM they may be 
enhanced; (iv) A large part of the parameter space so far survived will be soon covered by 
the expected XENON100(2012) sensitivity (especially for the NMSSM). 

Therefore, although the low energy SUSY can in general accommodate a SM-like Higgs 
boson near 125 GeV and enhance its di-photon signal rate at the LHC, not all models of low 
energy SUSY are equally competent if they are required to satisfy all current experimental 
constraints. From our present study and some other studies in the literature, we conclude: 



The fancy CMSSM/mSUGRA is hard to give a 125 GeV SM-like Higgs boson (l7|. 



The MSSM can give such a 125 GeV Higgs and can also enhance its di-photon signal 
rate at the LHC, which, however, will incur some fine-tuning. 

The nMSSM (the nearly minimal SUSY model) can give a 125 GeV SM-like Higgs, 
but severely suppress its di-photon signal rate at the LHC [6|. 

The NMSSM is so far the best model to accommodate such a 125 GeV Higgs; it can 
naturally (without fine-tuning) predict such a SM-like Higgs mass and readily enhance 
its di-photon signal rate at the LHC. At the same time, in a large part of its parameter 
space, this model can also enhance the signal rates pp ^ h ^ VV* {V = Z, W) at the 
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LHC and predict a large scattering rate of dark matter and nucleon at the XENONIOO. 
So the interplay of LHC and XENONIOO will soon allow for a good test of this model. 



Acknowledgement 

This work was supported in part by the National Natural Science Foundation of China 
(NNSFC) under grant Nos. 10821504, 11135003, 10775039, 11075045, by Specialized Re- 
search Fund for the Doctoral Program of Higher Education with grant No. 20104104110001, 
and by the Project of Knowledge Innovation Program (PKIP) of Chinese Academy of Sci- 
ences under grant No. KJCX2.YW.W10. 



[1] ATLAS Collaboration, larXiv: 1202. 1408) [hep-ex]; larXiv: 1202. 1414] [hep-ex]; ' arXiv: 1202. 14151 
[hep-ex] . 

[2] CMS Collaboration, larXiv: 1202.1416] [hep-ex]; larXiv: 1202. 1 487] [hep-ex]; iarXiv:1202.1488l 

[hep-ex]; larXiv: 1202. 1489] [hep-ex]; larXiv:1202.1997] [hep-ex]; larXiv:1202.3478] [hep-ex]; 

larXiv: 1202.36T71 [hen-ex]: laSCiv: 1202.4083] [hen-ex]: ' aSiv: 1202.4195] [hep-ex]. 
[3] D. Carmi, A. Falkowski, E. Kuflik and T. Volanskv. larXiv: 1202. 3 144] [hep-ph]. 
[4] J. R. Espinosa, C. Grojean, M. Muhlleitner and M. Trott. [arXrv:1202.3697] [hep-ph]: A. Aza- 

tov, R. Contino and J. Galloway, larXiv: 1202.34151 [hep-ph]. 
[5] S. Moretti and S. Munir, Eur. Phys. J. C 47, 791 (2006); K. Hsieh and C. P. Yuan, Phys. Rev. 

D 78, 053006 (2008); I. Low and S. Shalgar, JHEP 0904, 091 (2009); U. Ellwanger, Phys. 

Lett. B 698, 293 (2011). 
[6] J. Cao, Z. Heng, T. Liu, J. M. Yang, Phys. Lett. B 703, 462 (2011). 
[7] U. Ellwanger, larXiv:1112.3548i 

[8] H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985); J. F. Gunion and H. E. Haber, 

Nucl. Phys. B 272, 1 (1986) [Erratum-ibid. B 402, 567 (1993)]. 
[9] A. Djouadi, Phys. Rept. 459, 1 (2008) (arXiv:hep-ph/0503173l . 
[10] M. S. Carena, J. R. Espinosa, M. Quiros and C. E. M. Wagner, Phys. Lett. B 355, 209 (1995) 

|arXiv:h"ep^h/95043l6] . 

[11] P. H. Chankowski, J. R. Ellis and S. Pokorski, Phys. Lett. B 423, 327 (1998) |h ep-ph/9712234] ; 



19 



R. Barbieri and A. Strumia, Phys. Lett. B 433, 63 (1998) |h ep-ph/9801353] ; G. L. Kane and 
S. F. King, Phys. Lett. B 451, 113 (1999) |hep-ph/9810374]. 

[12] U. Ellwanger, C. Hugonie and A. M. Teixeira, Phys. Kept. 496, 1 (2010); M. Maniatis, Int. 
J. Mod. Phys. A25 (2010) 3505 |arXiv:0906.0777 [hep-ph]]. 

[13] For phenomenological studies, see, e.g., J. R. Elhs et al, Phys. Rev. D 39, 844 (1989); 
M. Drees, Int. J. Mod. Phys. A4, 3635 (1989); S. F. King, P. L. White, Phys. Rev. D 52, 4183 
(1995); B. Ananthanarayan, P.N. Pandita, Phys. Lett. B 353, 70 (1995); B. A. Dobrescu, 
K. T. Matchev, JHEP 0009, 031 (2000); R. Dermisek, J. F. Gunion, Phys. Rev. Lett. 95, 
041801 (2005); G. Hiller, Phys. Rev. D 70, 034018 (2004); F. Domingo, U. EUwanger, JHEP 
0712, 090 (2007); Z. Heng et al, Phys. Rev. D 77, 095012 (2008); R. N. Hodgkinson, A. 
Pilaftsis, Phys. Rev. D 76, 015007 (2007); W. Wang et al, Phys. Lett. B 680, 167 (2009); 
J. M. Yang, Int. J. Mod. Phys. D 20, 1383 (2011) j arXiv:1102.4942] [hep-ph]]: U. Ellwanger 
and C. Hugonie, Mod. Phys. Lett. A 22, 1581 (2007) [a rXiv:hep-ph/0612133 |; Eur. Phys. J. 
C 25, 297 (2002) [arXiv:hep-ph/9909260] ; U. Ellwanger, Eur. Phys. J. C 71, 1782 (2011) 
|arXiv:1108.0"T57] [hep-ph]]. 

[14] S. F. King, M. Muhlleitner and R. Nevzorov. [arXiv:1201.267T] [hep-ph]. 

[15] Z. Kang, J. Li and T. Li. larXiv:1201. 53051 [hep-ph]. 

[16] M. Bastero-Gil, C. Hugonie, S. F. King, D. P. Roy and S. Vempati, Phys. Lett. B 489, 359 
(2000) | hep-ph/0006"T98l ; A. Delgado, C. Kolda, J. P. Olson and A. de la Puente, Phys. Rev. 
Lett. 105, 091802 (2010) [arXiyilOOSJlHl] [hep-ph]]; U. Ellwanger, G. Espitalier-Noel and 
C. Hugonie, JHEP 1109 (2011) 105 [arXiv: 1107.24721 [hep-ph]]: G. G. Ross and K. Schmidt- 
Hoberg, arXiv:1108. 1284 1 [hep-ph]. 

[17] S. Heinemeyer, O. Stal and G. Weiglein, larXiv:1112.3026l A. Arbey et al, arXiv:1112.3028t 
L. J. Hah, D. Pinner and J. T. Ruder man . larXiv: 1112. 27031 P. Draper et aL. ' arXiv:1112.3068t 
A. Arbey, M. Battagha and F. Mahmoudi, larXiv:1112.3032| O. Buchmueller et al, 
larXiv:1112.3564 : M. Kadastik et al, larXiv:1112.3647l J. Cao, Z. Heng, D. Li and 
J. M. Yang, |arXiv:1112.4391l [hep-ph]. A. Arvanitaki and G. Vihadoro, arXiy:1112.4835J 
[hep-ph]; H. Baer, V. Barger, A. Mustafayev, larXiv:1112.30T7l I. Gogoladze, Q. Shafi 
and C. S. Un, larXiv:1112.2206 [hep-ph]; J. L. Feng, K. T. Matchev and D. Sanford, 
larXiv:1112:302r [hep-ph]; S. Akula, B. Altunkaynak, D. Feldman, P. Nath and G. Peim, 
larXiv:1112.3645 [hep-ph]; A. Bottino, N. Fornengo and S. Scopel, larXiv:1112. 56661 [hep- 



20 



ph]; J. F. Gunion, Y. Jiang and S. Kraml, larXiv:1201.0982l [hep-ph]; P. Fileviez Perez, 
larXiv:1201.150n [heD-Dh]: J. Ellis, M. K. Gaillard and D. V. Nanopoulos. larXiv:1201.6045] [heo- 
ph]; N. Karagiannakis, G. Lazarides and C. Pallis. [arXiv:1201.21 11' [hep-ph]; L. Roszkowski, 
E. M. Sessolo and Y. -L. S. Tsai. larXiv:1202.1503l [hep-ph]: L. Aparicio, D. G. Cerdeno and 
L. E. Ibanez. [a?Xiv:1202.0822l [hep-ph]: C. -F. Chang, K. Cheung, Y. -C. Lin and T. -C. Yuan, 
E-Xry:r202:005T ' [hep-ph]; K. A. Olive, arXiv: 1202 .23241 [hep-ph]; J. Ehis and K. A. Olive, 
larXiv:1202.3 262| [hep-ph]; H. Baer, V. Barger and A. Mustafayev, larXiv: 1202.40381 [hep-ph] ; 
D. Ghosh, M. Guchait and D. Sengupta. larXivri202.4937 [hep-ph]: N. Desai, B. Mukhopad- 
hyaya and S. Niyogi, larXiv:1202.5T90 ^ [hep-ph]. 
[18] M. Carena, S. Gori, N. R. Shah and C. E. M. Wagner. [a?Xiv:1112.3336] [hep-ph], 
[19] D. J. . Miher, R. Nevzorov and P. M. Zerwas, Nucl. Phys. B 681, 3 (2004) 
|arXiv:hep-ph /0304049] . 

[20] U. Ellwanger, J. F. Gunion and C. Hugonie, JHEP 0502, 066 (2005); U. Ellwanger and 

C. Hugonie, Comput. Phys. Commun. 175, 290 (2006). 
[21] G. Degrassi et al, Eur. Phys. J. C 28 (2003) 133; S. Heinemeyer, W. Hollik and G. Weiglein, 

Eur. Phys. J. C 9 (1999) 343; S. Heinemeyer, W. Hollik and G. Weiglein, Comput. Phys. 

Commun. 124 (2000) 76; M. Frank et al, JHEP 0702 (2007) 047; 
[22] G. Altarelh and R. Barbieri, Phys. Lett. B 253, 161 (1991); M. E. Peskin, T. Takeuchi, Phys. 

Rev. D 46, 381 (1992). 
[23] J. Cao and J. M. Yang, JHEP 0812, 006 (2008). 
[24] M. Davier et al, Eur. Phys. J. C 66, 1 (2010). 

[25] J. Dunkley et al. [WMAP CoUaboration], Astrophys. J. Suppl. 180, 306 (2009). 

[26] E. Aprile et al. [XENONIOO Collaboration], Phys. Rev. Lett. 107, 131302 (2011). 

[27] H. Ohki et al, Phys. Rev. D 78, 054502 (2008); D. Toussaint and W. Freeman, Phys. Rev. 

Lett. 103, 122002 (2009); J. Giedt, A. W. Thomas and R. D. Young, Phys. Rev. Lett. 103, 

201802 (2009). 

[28] N. Desai and B. Mukhopadhyaya, larXiv:lllL2830] [hep-ph]; X. J. Bi, Q. S. Yan and 
P. F. Yin, Phys. Rev. D 85, 035005 (2012) [a rXiv:llll. 22501 [hep-ph]]; B. He, T. Li, Q. 
Shafi, larXiv: 11 12.44611 [hep-ph]. 

[29] M. S. Carena, D. Garcia, U. Nierste and C. E. M. Wagner, Nucl. Phys. B 577, 88 (2000). 

[30] J. Cao and J. M. Yang, Phys. Rev. D 78, 115001 (2008). 



21 



[31] J. Cao, H. E. Logan and J. M. Yang, Phys. Rev. D 79, 091701 (2009). 

[32] G. Aad et al. [The ATLAS Collaboration], |ax2£ivi09QL0512J [hep-ex]. 

[33] J. Cao, et al, Phys. Rev. D 82, 051701 (2010). 

[34] J. Cao, et al, JHEP 1007, 044 (2010). 



22 



